The reproduction of the dynamics of a quantum system by an ensemble of classical 
particles beyond de Broglie— Bohmian mechanics 
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It is shown that for any given quantum system evolving unitarily with the Hamiltonian, 
H — p 2 /(2m) + U(q), [bold letters denote D-dimensional (D ^ 3) vectors] and with a suf- 
ficiently smooth potential t/(q), there exits a classical ensemble with the Hamilton function, 
^(P,q) = P 2 /(2m) + (7 (oo) (q), where the potential (7 (oo) (q) coincides with U(q) for almost all q 
(i.e., LA°°' can be different from U only on a measure zero set), such that the square modulus of the 
wave function in the coordinate (momentum) representation approximately equals the coordinate 
(momentum) distribution of the classical ensemble within an arbitrary given accuracy. Furthermore, 
the trajectories of this classical ensemble, generally speaking, need not coincide with the trajectories 
obtained from de Broglie-Bohmian mechanics. Consequences of this result are discussed. 

PACS numbers: 03.65.Ta, 03.65.Sq, 45.20.Jj, 02.30.Cj, 02.10.De, 02.30.Zz, 02.40.Pc 



I. INTRODUCTION 

It is well known that results of quantum mechanics can- 
not be reproduced by a local and/or context-independent 
hidden variable theory 111, Celebrated de Broglie- 
Bohmian mechanics [j-Q; being a non-local and contex- 
tual hidden variables theory, gives a consistent way of 
constructing an ensemble of classical trajectories which 
reproduces exactly quantum results. This is achieved by 
adding an extra term called the quantum potential, which 
depends on the wave function, into the Hamiltonian of 
the system. 

Is there any other systematic way (i.e., explicitly non 
de Broglie-Bohmian approach) of constructing an ensem- 
ble of classical trajectories that reproduces the quantum 
probability distribution? The answer to this question is 
affirmative owing to the recent method developed in Ref. 
that allows one to construct classical trajectories from 
quantum mechanics using solely the time evolving prob- 
ability density (not the full wave function) and without 
assuming or solving any equations of motion. Further- 
more, it follows from the construction that these trajec- 
tories need not coincide with de Broglie-Bohmian trajec- 
tories. The method in Ref. [7] employes the geometrical 
construction of centroidal Voronoi tessellations. 

In the current paper, we go further and present an al- 
ternative, though complementary, point of view on non 
de Broglie-Bohmian trajectories. Based on very mild 
topological assumptions on the configuration space, we 
show non-constructively (see Theorem [5]) that for any 
given quantum system evolving unitarily with the Hamil- 
tonian 



H = p 2 /(2m) + [/(q), 



(1) 



(bold letters denote /^-dimensional vectors throughout 
the paper) where the potential f(q) is a sufficiently 



smooth function, one can construct a classical ensemble 
with the Hamilton function 



Jf(p,q)=p 2 /(2m) + C/H( q ) 



(2) 



where the potential J7^°°^(q) coincides with U(q) for al- 
most all q (i.e., L/(°°) can be different from U only on a 
measure zero set), such that the square modulus of the 
wave function in the coordinate (momentum) represen- 
tation approximately equals the coordinate (momentum) 
distribution of the classical ensemble within an arbitrary 
given accuracy. Additionally, the quantum system with 
the Hamiltonian 



#(«>) =p7(2m) + m°° J (q 



(oo) 



(3) 
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is, physically speaking, the same as the original one with 
the Hamiltonian ([T]). 

Since our proof is non-constructive, the trajectories 
of the classical ensemble are only known to satisfy 
very basic geometrical (kinematical) assumptions; hence, 
these trajectories are not related to de Broglie-Bohmian 
ones in a general case. 

In a nutshell, the proof of Theorem [5] is based on the 
fact that on the one hand, a quantum system does not 
"feel" any measure zero set modifications of the potential 
due to the dynamical non-locality of quantum mechan- 
ics; on the other hand, the corresponding classical system 
"feels" the difference due to the dynamical local nature 
of classical mechanics. Here, the dynamical non-locality 
of quantum mechanics means that during unitary non- 
relativistic evolution, the value of the wave function in 
some point at a subsequent time moment depends on 
the values of the wave function and the potential in all 
the points at the current time moment. Correspondingly, 
the dynamical locality of classical mechanics implies that 
according to the Hamilton equations, the position of a 
classical particle at a subsequent time moment depends 
only on the values of the potential in some small neigh- 
bourhood of the current position of the particle (i.e., the 
subsequent position depends on the force that acts on the 
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particle in the current position). The measure zero set 
modifications of the potential U are done along curves, 
which shall be called "scratches," and the procedure of 
doing such modifications shall be named "scratching" the 
potential U. The main technique of the proof is to per- 
form scratching and locate a classical particle such that 
the scratches become trajectories of the particle. Then, 
the shapes of the scratches are chosen such that the prob- 
ability to find the classical particle in some region ap- 
proximates the analogous probability for the quantum 
particle. 

Strictly speaking, neither the Bell inequalities [l[ nor 
the Kochen-Specker theorem 0] applies to the classical 
ensemble constructed in Theorem [2] because Theorem [2] 
does not discuss the measurement process. A complete 
hidden variable theory can be obtained from Theorem [2] 
after embedding the measurement process. This should 
be done in future investigations; now however, we present 
intuitive arguments why such embedding should not con- 
vert Theorem [2] into a local or context-independent hid- 
den variable theory. As described in details in the proof 
of Theorem [2J scratching is performed depending on the 
actual values of the square modulus of the wave func- 
tions at different time moments and at different positions. 
Hence, the constructed classical ensemble is to violate 
the Bell inequalities by the same token as de Broglie- 
Bohmian mechanics does. The contextuality of Theorem 
[2] is to arise from the fact that the classical ensemble re- 
produces the square modulus of the wave functions only; 
it does not provide any information upon other observ- 
ables. 

The rest of the paper is organized as follows: The next 
section contains exact formulation and rigorous proof of 
the informally reformulated result. Discussions of physi- 
cal consequences of Theorem [2 such as the issue of the 
discrimination between classical and quantum dynamics 
based on a posteriori information, a new interpretation of 
quantum mechanics, and quantization and semi classical 
approximation, are presented in Sec. IIIH The conclu- 
sions are drawn in the last section. Finally, a result upon 
the Diophantine approximation that is employed in the 
proof of Theorem [2] is obtained in the Appendix. 



II. MAIN RESULTS 

Perhaps, a simplest result that can be obtained within 
our framework is Theorem [T] Despite the fact that Theo- 
rem[5]is the main result of the paper, we present Theorem 
[T] due to the following methodological motivation: The 
proofs of both the theorems are mainly based on same 
ideas, however, Theorem [T] allows for more lucid presen- 
tation of these ideas owing to a relative simplicity of its 
proof. 

Let A denote the Lebesgue measure in M. D throughout 
the paper. 

Theorem 1. Assume that D ^ 2 and 



1. B C MP , < A(£>) ^ oo, is an open convex set; 

2. B\ , . . . , B n is a closed cover of B such that B C 

l,...,n, A(5j n B k ) = 0, Vk^j; 

3. A quantum system with the Hamiltonian 

H(t)=p 2 /(2m) + U(q,t), (4) 

is given, where U (q, t) is a continuously differen- 
tiable function Vq £ B and a continuous function 
Vt £ [ijji/] (—oo < ti < tf < +oo) and 



[ \U{q,t)\d D q<oo, Vt £ [ti,t f ] 
J B 



(5) 



4- The wave function of the system, ^(q, t), evolves 
according to the Schrodinger equation \/t £ [ti,tf], 



[ |*(q,t)| 2 d D q=f, Vte[U,t f 
Jb 



Let 



Pk(t)= [ \^(<l,t)\ 2 d D q , 

Jb,. 



(6) 



(7) 



denote the probability to find a quantum particle in 
the region Bk at time moment t; 

Then, for an arbitrary natural number Q > n 2n , there 
exists an ensemble of N (0 < N ^ Q) classical particles 
with the Hamiltonian function 



^ (A) (P,q^)-P 2 /(2m) + [/W(q,t), 



(8) 



and a quantum system with the wave function (q, t) 
and the Hamiltonian 



m x \t)=p 2 /(2m) + U^{q,t), 



(9) 



where U^(q,t), VA > 0, is a continuously differentiate 
function Vq G B and a continuous function Vi £ 
such that the equalities 



lim U^(q,t) = tf(q,t), 



A— >OG 



lim <£ (A) (q,i) = *(q,i), 

A— too 



(10) 
(11) 



are valid Vt £ [ti, tf] and for almost all q £ B (they may 
be violated on a measure zero set). Furthermore, there 
exits the initial condition for the classical ensemble such 
that 

max { P k (Uj) -4 A) (*i,/) } < T7^. VA > 0,(12) 
fe=l,...,n I. J NQ^ 

where ir k ^ (t) denotes the probability to find a classical 
particle in the region B^ at time moment t, and the equal- 
ity 



5>< A) (t) = l, VA>0 



(13) 



fc=i 
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is valid for almost all t G [^i/] 1 , i.e., the motion of the 
classical ensemble is confined in the region B. 

Proof. The main motif of the current paper is the em- 
ployment of the procedure of "scratching" the potential. 
We call (q, t) the scratched version of the potential 
J7(q,t)if 



and Eq. f| 12[) is valid if we set 



N 

t/( A )(q,t) = (7(q,t)n [l 



-A/*(q) 



k=l 



(14) 



where /fc(q) is continuously differentiable, and the equa- 
tion /fc(q) = defines a differentiable curve, which is 
called the fc th scratch, and moreover, we assume that 
dfk/dqj — along the k th scratch. One readily notices 
that 



U {X \q,t) 





U(q,t)+0(X-°°) 



Vf=x[/*(q) = o] tl5) 

otherwise, 



where \J is the logical disjunction and the symbol 
O (A~°°) denotes a term that decays faster than any 
power of 1/ A as A — > oo. Equation (fT5j) can be rephrased 
as the potential along the scratches vanishes. 

One of most important property of the procedure of 
scratching is that 



017(A) 
dq.j 





dU/dqj + 0{X~°°) 



V£=i [/*(q) = o], 

otherwise. 



(16) 



This equation implies that a classical particle experiences 
no force along the scratches. Using this idea, we force N 
classical particles to move along straight lines by choosing 

/fc(q) = (q-q*(*i)-{q*(t/)-q*(ti)} (17) 

x {[qJi-MtiM/toCt/Jli-MtOh}) 2 , 

where [a] j denotes the j th component of a vector a. Since 
/ fe (q) = <=> q = q fe fo) + [q fe (t/) - q*^)] & (18) 
£ = {[q]i - [qfc(^)]i}/{[q*(t/)]i " [q*(«i)]i> , 

and also £ (<!*(£»)) = 0, £ (q/.(i/)) = 1, we indeed con- 
clude that equation /fe(q) = defines the line that 
goes through the points qfc(ti) and (\f.(tf); additionally, 
dfk/dqj = on that line. 

According to Eq. ©, YT k =i p k{U) 

Y^k=i-Pk{tf) = 1; hence, the sequence of real 
numbers P\{U), ■ ■ ■ , P n {U), Pl(tf), ■ ■ -,P n (tj) 

obeys the conditions of the Lemma presented 
in the Appendix. Therefore, we conclude that 



VQ G N, 



> n 



2 1 1 



there exists (positive) integers 



,JV„(ti),JVi (*/),. 



fc=i 



E 

fc=i 



, N n (tf), such that 



1 There may be time moments when the classical particle is on 
the boundary of one of the regions Bj. . 



,(A) 



(Uj) = N k (Uj)/N. 



(20) 



According to Eqs. ([T5]l and (f^Uj) . it is justifiable to inter- 
pret N as the total number of particles in the classical 
ensemble, and Nk{U) [Nk(tf)] as the number of classical 
particles located in the region B k at time moment U [tf]. 

Note that there are in fact two possible interpretations 
for the positive integer N. First, we could assume that 
we have N identical copies (i.e., N realizations) of the 
one-particle classical system with the Hamiltonian ©. 
(This is a usual construction in literature.) Obviously, 
the particle from one such a copy cannot interact with 
the particle from another copy. Second, the interpreta- 
tion utilized here as well as in Theorem [5] is that we as- 
sume to have one A^-particle classical system (i.e., "gas" 
of N particles), where these particles can interact only 
when they collide with each other. Since we will locate 
scratches and select the initial conditions to completely 
eliminate collisions, both these interpretations are indeed 
equivalent. Nevertheless, the latter one is more general. 

Let us arrange the motion of the N particles accord- 
ing to the interpretation. This, however, can be easily 
achieved by means of the scratched potential [Eq. (fl"4"l) ] 
with Eq. (fl7|) . where qk(ti.f) denote the initial and final 
positions of the fc th particle, correspondingly. Note that 
each particle must be placed in its own scratch. The final 
problem is to select points (tk{U,f) such that no collision 
between particles is possible. 

Let us construct the finite sets Ckitij), k = 1, . . . ,n, 
such that C k (ti,f) C B k \dB k , \C k {Uj)\ = N h (U,f), 

71 

(J C k {Uj) = {qi(tij),...,q N (tij)}, 
k—l 

and the following property is satisfied: 



Vx, y, z G C Vt G 



z^x + i(y-x), (21) 



where C = ULi C k(U) U C k {t f ) and |C| = 2N (here, 
\A\ denotes the number of elements of a finite set A). In 
other words, C is a set of all the initial and final positions 
for the classical particles; Ck(U) [Cfc(i/)] is a set of all 
initial [final] positions for particles in the region B k . All 
IN elements of the set C are distinct points, and no three 
points from C lie on a line. Such selected points Qk(tij) 
guarantee not only that no two particles share a scratch, 
but also that the collisions are prevented (nevertheless 
only partially). 

We shall prove by contradiction that the described 
above set C exists. The justification of the existence 
of C is based on the fact that A(Bk) ^ 0, k = 0, . . . , n. 
Let us fix integers k,l,m = l,...,n (k, I, and m may 
coincide) and introduce the notation A^ = Bk\ dBk, 
=Bi\ dB u and = B m \ dB m . The statement 

3x G A {1 \ y G B {1) Vz G C* (1) 3t G R z = x + t(y - x) 
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must be false because otherwise it implies that the set 
is a line, hence A (C^ 1 )) = 0, which contradicts as- 
sumption 2. Therefore, the following statement is true: 

Vx e A (1 \y e 3z g C (1) \/t e E z ^ x + t(y - x) 

Let us pick an arbitrary triple of such points and denote 
them by x' 1 - 1 , y*- 1 -*, and z^, respectively. Now construct 
the sets = \ {x^}, B^ = \ {y^}, and 
C (2) = C (i) \ | z (i)}. The statement 

3x G A (2) ,y e 5 (2) Vz 6 C (2) 3i G E z = x + i(y - x) 

is false because A(C (2) ) = A(C (1) ). One readily ob- 
serves that this iteration procedure can be continued an 
arbitrary number of times; hence, the sets Ck(Uj) exist 
for any N. 

Recall that qfc(ij,/) denote the initial and final posi- 
tions of the k th particle. If the initial momentum of the 
k th particle is chosen as 

p fc = m [qjt(tf) - qfe(ti)] /{tf - U), 

then, according to the Hamilton equations, its trajectory 
is 

q*(f) = qfe(*i) + {t- U) [qje(tf) - qk{k)} /{tf - U). 

Assume that the j th and I th particles collide, viz., 3t G 
suc h that qj(£) = (£) , i.e., 

t-tj = [<b(tj) - q;(^)]i 

t f - U [qi(t f ) - qj(t f ) +qj(U) - miU)^ 
= faj(U) -cu(U)] D 

[qj(*/)-qf(t/) + qj(tO-<«(ti)]£> 

However, because B is a Hausdorff space, we can always 
perturb, say [q :) (t i )] 1 , by some small value such that this 
chain of equalities is violated, but the sets Ck(t%j) pre- 
serves their properties. 

We note that Eq. {T5J follows from the fact that the 
trajectories are linear and the convexity of B (assumption 

The Schrodingcr equation in the momentum represen- 
tation for the scratched quantum system with the Hamil- 
tonian © reads 

i»^$W(p J t) = *W(p J t)p 2 /(2m) 

+ / d D k[7w(k-p,<)$( A )(lM), (23) 
Jr d 

where 

<& (A) (p,i) = J B {2 f h */2 tt (A) (q,*)e- ip " q/R , (24) 

^ )(P '* } = Sb (2^ C/(A)(q ' <)etP q/ ' 1 - (25) 



The integral over k in Eq. (|23f) manifests the dynamical 
non-locality of quantum mechanics. From Eq. (|15j) and 
the fact that the value of an integral is unchanged if the 
integrand is modified on a measure zero set, we conclude 
that 

UM{p,t) = U{ V ,t) + 0{\-°°). (26) 

Assumption guarantees the existence of the Fourier 
transform of U as well as U^ x \ Physically speaking, Eq. 
(|26p is a direct consequence of the dynamical non-locality. 

Equations (|10|) and (fTT|) follow from the fact that the 
Fourier transform is "insensitive" to measure zero set 
modifications. □ 

Theorem 2. Assume that D ^ 3 and 

1. B C R D , < A(B) ^ oo, is an open path connected 
set; 

2. B\ , . . . , zs a closed cover of B such that B C 
U; =1 i^ A(B) = £; =1 A(5,) 7 A(fl 3 -) ^ 0, j = 

A(B i nB k )=0, Vfc^j; 

5. ^4 quantum system with the Hamiltonian 

H = p 2 /(2m) + [/(q), (27) 

is given, where U (q) &emg a twice continuously dif- 
ferentiate function Vq G B, and 

E7(q)>0, VqeB, / |?7(q)| rf^q < oo; (28) 

^. TTie wave function of the system, ^(q, i), evolves 
according to the Schrddinger equation Vi G 
(—00 < U < tf < +oo) and 

f \9(q,t)\ 2 dPq=l, Vt€ [<,,«,]. (29) 
J B 

Let 

P k (t)= f |*(q,t)| 2 ^q, (30) 
JB k 

denote the probability to find a quantum particle in 
the region B k at time moment t; 

5. Bi,...,B n is a closed cover of M. D such that 
^ D = U"=i-B„ ; A (5,) ± 0, j = l,...,n, 
A (fij n Bkj = 0, Vfc ^ j, and 

p k (t)= f m P ,t)\ 2 d D P , (3i) 

JB k 

being the probability that the momentum of the 
quantum particle is in the region B k and <£>(p,£) 
is the wave function in the momentum representa- 
tion; 
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6. U = ti < < 2 < • ■ • < *k-i < 



t 



f> 



Then, for an arbitrary natural number Q > n 2Kn , there 
exists an ensemble of N (0 < N ^ Q) classical particles 
with the Hamiltonian function 



^ (A) (p,q) = p7(2™) + t/ (A) (q), 



(32) 



and a quantum system with the wave function \E , ( A )(q, t) 
and the Hamiltonian 



P 2 /(2m) + f/( A )(q), 



(33) 



where U^(q), VA > 0, is a twice continuously differen- 
tiate function Vq S B , such that the equalities 



lim U W (q) = U(q), 
lim ¥ x \q,t) = *(q,t), 



(34) 
(35) 



are valid Vt G [ti,tf] and for almost all q e i? (^gs. Jff^| ) 
and 15*5)) may 6e violated on a measure zero set). Fur- 
thermore, there exits the initial condition for the classical 
ensemble such that 



max 
fe=l,...,j 



max 
fe=i,...,) 



{\Pk(t 3 

{\h(tj 



Vj 



2A'ii 



(36) 
(37) 



where 71"^ (t) denotes the probability to find a clas- 
sical particle in the region Bk at time moment 



t 



r J»>, 



'(*) 



lim > 



(t) and 7r[ A ^ (t) denotes 
the probability that the momentum of the classi- 
cal particle is in the region Bk at time moment t 
^°°' > (t) = lim.\->.oo Trl^ (t) . Additionally, the equalities 



(38) 



are m/id /or almost all t £ [ti,tf] 2 , i.e., the motion of 
the classical ensemble is confined in the regions B. 

Proof. The main idea of the proof is, again, the scratching 
procedure. Assume that all the scratches given implicitly 
in the region B, i.e., the set of solutions of the following 
system of equations 



F 1 W (q) = 0,. 



• ,FiL(q) 



o, 



(39) 



where each function is continuously differentiable, defines 
the / th scratch (I = 1, . . . , N). We assume further that 



2 There may be time moments when the classical particle (the 
momentum of the classical particle) is on the boundary of one of 
the region B k (B k ). 



there exists integer j (1 ^ j ^ D) such that the Jacobi 
matrix 



(1) 



D-l 



d(qi,...,qj-i,q j+ i,...,q D ) 



is non singular on the I th scratch. We point out that it is 
a natural assumption because according to the implicit 
function theorem, if Ji(j) is non singular then Eq. (|39j) 
defines (locally) a curve. 
Defying auxiliary functions 



D-l 

/*(q) = £ [if(q) 

3=1 



(40) 



we introduce the scratched potential 

N 

C/ (A) (q) = C/(q)]j[l-e- A/!(q) 
1=1 



(41) 



One readily verifies from definition (|40| that /;(q) = 
and dfi (q) /9q = if and only if q lines on the I th scratch. 
Setting Fp\q) = 0, I = 1, . . . , N, let us scrutinize the 
Hessian matrix of the function /; (q) on the Z th scratch, 



H(fl)\i = 



( d 2 h 



\ dq r dq k 
= 2 (x« 



D dF (l) dp (l) 

^-J dq r dq k 



(42) 
(43) 



where x 



(0 



and 



Ji = 



(dF®/dq k ,...,dF$/dqk j 
d(F[ l \...,F^)/d( qi ,...,q D ) 



On the one hand, according to Eq. (|42p. H(fi)\ l is a 
Gramian matrix; thus, H(fi)\ l is positive semidefinite. 
On the other hand, Eq. (|4"3"|) manifests that the Hessian 
H(fi)\, is of D — 1 rank because J; is singular (since 
F$ = 0) and the following D - 1 minor of the Hessian: 



d 2 h 



r ^ 3 



is not singular (since Ji(j) is non singular). Hence, we 
conclude that the Hessian H(fi)\ l has one zero eigenvalue 
and D — 1 positive eigenvalues. 

We now show that the eigenvector corresponding to 
the zero eigenvalue is the tangent vector of the scratch. 
Let q = q^(s) be a parametric representation of the ^ th 
scratch. Then, 



F, 



f (q%)) 



= 0, 



■ F, 



(0 



D-l 



q«( S )) =0 
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and consequently, dF^ /ds = 0, . . . , dF^'^/ds = 0. 
Since 



inequalities (1361) and (|37|) hold if 



1 ffCft) 1 ' 1 



D-l 



dF. 



= 2 E ^ 



0, (44) 



and the fact that the vector dqyi/ds is proportional to 
the tangent vector of the Z th scratch, we have confirmed 
our statement. 

Equations (fT5|) and (|16l) are also valid in the case of 
the potential (l4"Tj) . In other words, the scratches are equi- 
librium positions. Therefore, to analyze dynamics of the 
classical system (f3"2j) , we can linearize the Hamilton equa- 
tions in a neighbourhood of the I th scratch (such a pro- 
cedure is justified by the theorem on page 100 of Ref. 

!)• 

The Hessian of the scratched potential (|4"Tj) on the / th 
scratch reads 



H 



N 



XU H(f, 



l)\i 



n 



).(45) 



To understand the qualitative behaviour of the classical 
system, one needs to analyze the eigenvalues of the Hes- 
sian ([15]) (see the theorem on page 104 of Ref. [8]). Re- 
calling the properties of H(fi)\ l and inequality ([25]) . we 
conclude that the Hessian (|45|) has D — 1 positive eigen- 
values and one zero eigenvalue, with the corresponding 
eigenvector being the tangent vector of the I th scratch, if 
no scratches intersect and A > 0. 

According to the theorem on page 76 of Ref. [8] , if a 
classical particle is placed in a scratch, then in the limit 
A — > oo, the motion of the particle is constrained to the 
scratch and its velocity can be only collinear to the tan- 
gent vector of the scratch. We are now in the position to 
formulate the fundamental property of the scratches: the 
procedure of scratching a potential of some unconstrained 
classical mechanical system (implicitly) introduces holo- 
nomic constrains in the limit A — > oo. However, these 
constrains can be "felt" by a classical particle only if it 
is placed in one of the scratches, thus, the name implicit 
holonomic constrains. 

The same conclusion can be reached intuitively as fol- 
lows: If a classical particle, being initially in a scratch, 
does not move along the scratch, then according to Eq. 
(|4"5)) . the particle experiences the force of the magnitude 
proportional to A that pushes the particle back to the 
scratch. Hence, the motion of the particle evidently be- 
comes constrained in the limit A — > oo. 

The sequence of real numbers 

Pi(ii), . . . , P n (h), . . .,Pi{t K ), . . . , P n (t K ), 
Pi(ii),...,P„(ii),...,Pi(t K ),...,P n (t K ) obeys the 
conditions of the Lemma from the Appendix because 

Therefore, VQ £ N, Q > n 2Kn , there exists (positive) 
integers N,Ni(ti), . . . ,N n (h), . . .,Ni{t K ), ■ ■ .,N n (t K ), 
N 1 (t 1 ),...,N n (ti),...,N 1 (t K ),...,N n itK) such that 



r (°°) 



(t s ) = N k (tj)/N, 

3 = 1, 



,K, 



k = l,...,n. (46) 



Similarly to Theorem[TJ we interpret N as the total num- 
ber of classical particles in the ensemble, Nk(tj) - the 
number of classical particles in the region B k at time 
moment tj, and N k (tj) - the number of classical par- 
ticles with momenta in the region B k at time moment 
tj. 

For the given natural numbers Nk(tj), we can 
construct simple (i.e., non self-intersecting) paths 
qW(s),...,qW(s) : [0,1] -> B, such that there exist 
real numbers, — si < S2 < • • • < 8K = lj and 

|{qW( Sj ), ■ ■ ■ , q (A0 (^)} n ^ \ dB 4 = N ^)M7) 

where k = 1, . . . , n and j — 1, . . . , K. Moreover, we put 
an additional requirement that no such two paths can 
intersect. We shall demonstrate that such a construction 
exists. 

Due to the path connectedness of B, we can pick any 
finite set of distinct points in B and connect them by a 
path. Now we demonstrate that since B is an open set, 
a path connecting a finite set of distinct points can be 
chosen to be a simple one. Let us select any path that 
goes thought these points and denote it by q : [0,1] — > 
B. Assume that the path self-intersects in the point x. 
Since B is open and Hausdorff, we can select an open ball 
centred at the point of self- intersection, O e (x) = {y £ 
R D | ||y - x|| < e}, such that O e (x) C B and the radius 
of the ball is chosen such that the following condition 
satisfies 

{q(s) |0 < a *S l}f|90 e (x) = {q(si),q(s 2 ),q(s 3 ),q(s4)h 

where s% < S2 < S3 < S4 and q(si) 7^ q(sa) 7^ 
q( s 3) 7^ q( s 4)- In other words, the open ball contains 
two segments G\ = {q(s) | s\ < s < s 2 } C O e (x) and 
G*2 = {q(s) I s 3 < s < S4} C O s (x) which intersect 
Gi n G2 = {x}. We shall use the convexity of O e (x) 
and replace by linear segments that do not inter- 
sect. Let L(a, b) = {a + t(h — a) | ^ t ^ 1} denote a 
linear segment connecting points a and b. The following 
statement is true 

Vbi.ba.ei.ea £ dO £ (x) (bi ± b 2 + e x + e 2 ) 3c £ O e (x) 
[L(bi, x) U L(x, ei)] f) [L(b 2 , c) U L(c, e 2 )] = 0, 

because its negation contradicts the fact that O e (x) is 
an open ball (note that the assumption D ^ 3 is cru- 
cial). Substituting Gi j2 by L[q(si),x] UL[x, q(s 2 )] and 
L[q(s3), c] UL[c, q(c4)], correspondingly, we have demon- 
strated that there exists a simple path that connects a 
finite number of distinct points. 

Now we shall demonstrate that B' = B \ p(a, b) is 
an open path connected set, where p(a, b) = {q(s) | ^ 
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s ^ 1} C B denotes a simple path connecting a, b G B. 
Construct the set W E = {y G R D ||y - q(s)|| < e, < 
s ^ 1}. Since B is open, we fix e > such that dW £ C B. 
Let us choose arbitrary points c, d G B' . There are two 
possibilities - either p(c,d) C B 1 or p(c,d) <£_ B'. The 
latter case means that p{c, d) C B and p(c, d) DdW £ > ^ 
for some e' (0 < e' ^ e). However, all points in the 
set p(c, d) n dW e i can be path connected with all the 
paths lying in dW e < because dW e < is a path connected 
set. Finally, since dW e > C B' , we have demonstrated 
that B' is an open path connected set. By the same 
token, Vf , g e B', B" = B' \ p(f , g) (where p(f , g) C B ) 
is an open path connected set as well. In other words, we 
can subtract any finite number of simple non-intersecting 
paths from B without affecting its path connectedness. 
Therefore, qW(s), . . . , q( N \s) exists with the prescribed 
properties. 

The paths q^(s) are continuous curves by defini- 
tion. Thus, due to the Weierstrass theorem, they 
can be approximated polynomially within an arbi- 
trary given accuracy [and preserving property (|47|) 



Let qW(s), . . . , q( N \s) denote these polynomial curves. 
Having alter the I th curve in an arbitrary small neigh- 
bourhood of the point q^ l '{sj), Vj = 1, . . . ,K and V7 = 
1, . . . , N, such that 



dqW(«j)/d* 



dq(0( Sj )/d' 



A k = 



j\ f| A * = N k(tjl (48) 

" ) 1 = 1, ...,N 

x ^ and x £ B k \ dB k \ (49) 



(viz., the direction of the tangent vector of q^(s) at s = 
Sj coincides with the direction of one of the momentum 

vectors from the set Bk), we can use these new curves 
as the scratches, qW(s), . . . , qW(s), to enforce desired 
dynamics of the classical ensemble. 

Now let us slightly generalize the procedure of scratch- 
ing the potential U(q) by performing the following sub- 
stitution 



N 

E/W(q) = t/(q)-Q[i 



JV 



1=1 



i=i 



-A/i(q) 



(50) 



The fundamental property of the scratches, derived orig- 
inally for the case of Eq. (l4"Tj) . is valid in the case of 
the modified scratching [Eq. (|50l) ] if the potential Vi is 
a function of the curvilinear coordinate that is in the di- 
rection of the Z th scratch (see the theorem on page 76 of 
Ref. |8[). Here, the requirement that the scratches must 
be simple curves that do not intersect leads to an impor- 
tant conclusion that all the functions Vi can be defined 
Vq G B (evidently that the function Vi cannot be defined 
at points where two scratches intersect and at points were 
a scratch self- intersects) . 

Now we employ the formalism of the Lagrange equa- 
tions of the second kind (see, e.g., Ref. Q). Assuming 
that the parametric definition of the Z th scratch reads 



q = qW (s^'), it is natural to use as the generalized 
coordinate by interpreting this parametric form of the 
scratch as a holonomic constrain. The Lagrange equa- 
tions of the second kind for the generalized coordinate 
is as follows 



ni- 
di 



/dqW\ 2 dsW 
V da(9 J dt 



dVi 



l=l,...,N. (51) 



The final step of the proof is to construct the potentials 
Vi such that the interpretations of (tj ) and 7r^°°' (tj ) 
given by Eq. (|46p are correct. If we require that 



,(') 



[tj) 



3 ' 



(52) 



then according to Eq. (j47|). dynamics of the classical 
ensemble indeed realizes the interpretation of 7r[,°°'(^j). 
Recall that the momentum of the I th classical particle is 



p(%) 



dq« 

l : 

dt 



dq« ds® 



'dsO dt y j 



Since we have already chosen suitable directions of the 
momenta in Eq. (|48p . we are only left to find proper 
magnitudes of the momenta. Vx G 3c G K cx G 
Bk \ dBk; whence, there exit G R such that 



mci 



f| (B k \ dB k 

1=1,. ...N V 

j = l,...,K, fc = l,. 



Nk(tj)(53) 



If ds^ l \tj)/dt = cij, then the interpretation of {tj) 
is realized. 

Therefore, in order to finalize the proof, we need to 
find the potentials Vi by knowing that s"'(tj) — Sj and 
ds">{tj)/dt = c/j-. These type of problems are known as 
inverse problems for ordinary differential equations, and 
as demonstrated in Ref. (see theorem 4 of Ref. [Tcj ) , 
the potentials Vi can be found such that s^{tj) ~ Sj and 
ds"'{tj)/dt ~ Q.j within an arbitrary given accuracy. 

Conditions (f38|) follow from the construction of the tra- 
jectories of the classical particles. Equations and 
(|35|) are proven in the same manner as Eqs. (|10p and 
dTTJ) from Theorem [TJ □ 



III. PHYSICAL INTERPRETATIONS OF THE 
RESULT 

Main physical consequences of Theorem [2] are the fol- 
lowing: 



A. Is the system in the black box quantum or 
classical? 

Imagine that some experimentalist is given a black box, 
and he or she needs to determined whether a quantum 
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system or an ensemble of classical particles is inside the 
box. The Hamiltonian (i.e., the Hamiltonian operator 
in the quantum case and the corresponding Hamilto- 
nian function in the classical case) of the system inside 
is assumed to be known a priori. The experimentalist 
is only allowed to measure the coordinate and/or mo- 
mentum distributions at arbitrary time moments. Then, 
Theorem [2] negatively answers this question, i.e., having 
only coordinate and/or momentum distributions for some 
Hamiltonian system, one cannot conclude whether they 
are obtained as the squared modulus of a wave function 
or as distribution functions of an ensemble of classical 
particles. 

Let us explain how such a conclusion is reached from 
Theorem [21 The set B in the theorem in fact represents 
the given black box; thus, we assume that the topology 
of the black box obeys assumption 1 of Theorem [2] (in- 
deed, this assumption is valid for a majority of realistic 
experimental setups). It is crucial that no-matter how 
precise experimental equipment is, there is aways a finite 
accuracy in each and every experimental measurement. 
Additionally, no probability distribution of a continuous 
variable (such as coordinates and momenta) can be mea- 
sured, strictly speaking. What one measures in such a 
situation is the probability that a value of the variable 
lies in some small (but finite) region. In terms of The- 
orem [21 this means that neither ^(q, t)\ nor |$(p, t)\ 
is measurable; however, both Pk(t) [Eq. ([30]) ] and Pk{t) 
[Eq. ([31]) ] are measurable in a real-life experiment. At 
this moment, it is convenient to informally reformulate 
Theorem [2] from the point of view of the experimental- 
ist: Say that he or she first assumes that the system in 
the black box is quantum, then the collected data ought 
to be interpreted as Pk(t) and Pk(t), correspondingly. 
(Note that any real device is capable of taking only a 
finite number of "shots," so t is a discrete parameter.) 
However, the theorem says that the original potential 
energy can be slightly modified such that no real (i.e., 
of a finite precision) experiment can detect a difference 
[Eqs. ([M)l and ([331) ] between the original and altered 
potential. Moreover, one can construct an ensemble of 
classical particles (with the Hamiltonian function of the 
natural form where the altered potential energy used as 
a potential term [Eq. (|32l) ]) which reproduces the mea- 
sured data Pk(t) and Pk(t) within an arbitrary desired 
precision [Eqs. ([3"o]) and ([57])]. Thus, the negative an- 
swer. 



The question whether a given dynamical process being 
quantum or classical is of active interest in the field of 
quantum control (see, e.g., Ref. [Til] ) - Reiterating our 
results, we conclude that almost any given quantum evo- 
lution can be reproduced by means of a classical ensemble 
with a scratched potential; even the interference in the 
double-slit experiment can be obtained in this way. 



B. Theorem [2] as a new interpretation of quantum 
mechanics 

De Broglie-Bohmian mechanics gives a consistent way 
of constructing a classical ensemble which reproduces ex- 
actly quantum results. This is achieved by adding an 
extra term, which depends solely on the wave function, 
into the potential energy of the system. Theorem [2] gives 
an alternative way of construction a classical ensemble 
that reproduces quantum results, however, not exactly 
but within an arbitrary given accuracy. Contrary to the 
scratching procedure introduced in the current paper, 
the de Broglie-Bohmian modifications of the potential 
are quite "noticeable," i.e., they are done not on a mea- 
sure zero set; furthermore, the obtained classical ensem- 
ble is unconstrained. As it was discussed in the proof of 
Theorem [2] (see the fundamental property of scratches), 
the scratches are in fact holonomic constrains. How- 
ever, they are not usual constrains (we used the name 
"implicit holonomic constrains" ) since a classical parti- 
cle must have peculiar initial conditions in order to be 
able to experience their presence. 

In summary, an ensemble of classical particles can sim- 
ulate any quantum system if the potential energy of the 
system is properly adjusted (according to Bohmian me- 
chanics), or if (implicit holonomic) constrains are as- 
signed to the system (according to Theorem [2]). Even 
though the scratching procedure is dual to de Broglie- 
Bohmian mechanics, they both share a common feature 
- they are nonlocal hidden variable theories. The hidden 
variables in our approach are the scratches that cannot 
be experimentally detected due to their "measure zero 
set" nature. 



C. Quantization and quasi classical approximation 

Before presenting the connection between Theorem [2] 
and the procedures of quantization and quasi classical 
approximation, we recall the definition of the Hilbert 
space Jz? 2 (-B). The axiom of quantum mechanics 
postulates that elements in ££ 2 (B) represent quantum 
states. First, we construct the vector space L 2 (B) — 
{/ : B — y C | / is measurable and |/| 2 integrable} , 
where measuring and integration is done with respect 
to the Lebesgue measure. The sesquilinear form 
(fid) — Jb f( c l)9( c 0^(^ c 1S n °t an i nner product on 
L 2 (B) because the equality (/, /) = implies that 
/ vanishes almost everywhere (a. e.), but does not 
necessary mean that / = 0. Having introduced the 
space, jV = {/ 6 L 2 (B) | / = a. e. on B}, we define 
Jz? 2 (B) — L 2 (B)/,yY as a factor space. One can demon- 
strate that ([/],[<?]) = J B /(q)^A(dq), / £ [/] and 
g G [g], is indeed an inner product in ^f 2 (B). Here, [/] 
and [g] are equivalence classes of functions in L 2 (B). 

In the light of the presented above definition, Eq. ([35j) 
implies that ty^ 00 *) and ^ correspond to the same physical 
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state. By the same token, Eq. (|34j) means that, as far 
as quantum mechanics is concerned, the potentials [T^ 00 -* 
and U correspond to the same quantum system; however, 
they correspond to two qualitatively different classical 
systems. 

Hence, generally speaking, quantization maps a classi- 
cal system to an equivalence class of quantum systems 
modulo scratching. Since the procedure of quantiza- 
tion and quasi classical approximation are closely related 
[l2l [l3| , the "dual" statement is also true: quasi classical 
approximation maps a quantum system to an equivalence 
class of a classical systems modulo scratching. 

We reiterate that scratching may tremendously change 
the behaviour of a classical system. For example, con- 
sider a system that contains a potential barrier such that 
a classical particle with some fixed energy cannot go over 
the barrier. The quantum counterpart of the classical 
particle of course can "go through" the barrier. To em- 
phasize the conceptual difference in the behaviours of the 
classical and quantum systems, we use the term "quan- 
tum tunnelling." Nonetheless, one can readily scratch 
the barrier such that the classical particle can also "go 
through" the barrier, but at the same time, this alter- 
ation of the potential has no effect on the behaviour of 
the quantum counterpart. 

Physically speaking, the fact that scratching does not 
affect the behaviour of a quantum system is due to the 
dynamical non-local character of quantum mechanics. 
Correspondingly, the dynamical local character of classi- 
cal mechanics is responsible for the effect that scratching 
has on a classical system. 



IV. CONCLUDING REMARKS 

In Theorem [21 we have demonstrated that a classical 
ensemble can be constructed that reproduces within an 
arbitrary accuracy the coordinate and momentum prob- 
ability distributions of a given quantum system, and the 
potential energy terms in the Hamiltonians of the clas- 
sical and quantum systems almost coincide. The con- 
sequences of the existence of such a construction have 
been discussed. It is important that the trajectories of 
the classical ensemble may differ from the trajectories 
obtained from dc Broglic-Bohmian mechanics. 

Theorem [5] complements results of recent paper Q. In 
fact, Theorem [2] can be proven constructively employing 
the geometrical construction of centroidal Voronoi tessel- 
lations in the same manner as it was utilized in Ref. Q. 
In such a case, the regions Bj from assumption 2 of The- 
orem [2] should be substituted by the Voronoi cells, Cj, 
see Eq. (2) of Ref. Q; the points qW(tj) should be sub- 
stituted by the centroids of the corresponding Voronoi 
cells, x/ = Xi(tj), see Eq. (3) of Ref. [7j|. Then, the tra- 
jectories constructed by means of the centroidal Voronoi 
tessellations can be used as scratches. However, note that 
the current proof of Theorem [5] is far more general than 
the summarized approach because it solely relies on very 



basic topological assumptions. 

We want to emphasize the importance of the topology 
of B in Theorem [51 i.e., the region of the configuration 
space where the quantum system is confined. The re- 
gion B represents the geometry of a measuring appara- 
tus (viz., the black box in the terminology of Sec. IIII Al) . 
Had the region B been a disconnected topological space, 
the idea of scratching, generally speaking, could not have 
been utilized with the requirement that the classical en- 
semble should also be confined in B. 3 

The property that the classical ensemble reproduces 
the momentum distribution is not a primary one since 
the measurement of any observable can be reduced to 
measuring coordinates. As a matter of fact, this prop- 
erty can be completely removed from the formulation of 
Theorem [3 without affecting the rest of the statement; 
in this case, steps (|48|) and (l53l) need to be removed from 
the poof, and the upper bound on Q can be relaxed to 



> n 



Kr. 



Theorem[5] allows for many further generalizations. For 
example, it would be also interesting to apply Valentini's 
analysis 0, [H[ to Theorem [2] in future papers. 
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Appendix A: A result regarding the Diophantine 
approximation of a sequence of real numbers with 
linear constrains 



Lemma. Vaf , . . . , an\ ■ ■ ■ , a[ K \ . . . , a„ 6 
that 



such 



5>} r) = aW/bW, A^lB^eZ, r = l,...,K, 
3=1 



VQ 



N, Q (nmax r {|B( r )|})" K , 



3q, , . . . , an' , ■ ■ ■ , oi" ; , ■ ■ ■ , an" ; £ Z, such that 

\\a^-a^/q}<-^, 0<g<Q,(AI) 



.(!) 



max 

3=1,. ..,n 
r=l,...,K 



qQ- 

n 

Y / a\ r) /q = A^/B^, r=l,...,K.(A2) 

3=1 



3 There is, nevertheless, an exception. If B is disconnected, 
i.e., B = B(l) U B(2) and B(l) n B{2) = 0, then one might 
expect that Theorem [2] should be valid without any changes 
if / s{1) |vK(q,t)| 2 ^q = C and / B(2) |*(q,t)| 2 d^q = 1 - C, 
Vt 6 [ti,tf]. In this case, we can apply the scratching proce- 
dure in each subregion B(l) and B(2) separately since there is 
no "exchange" of the probability between them. 
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Proof. According to the Dirichlet approximation theorem Assuming that Q > (n max r { | B M | } ) nK , we obtain 
(see also Ref. [13]), Vaf ) , . . . , a4 K) e K, VQ G N, 



Q > 1 , 3g, aj , . . . , a« £ Z, such that 



max < 

i=l,...,n L 



(r) (r) 



Therefore, Eq. (|A1[) is proven. To derive Eq. (|A2I) . 
consider 



Ar) 



BW 



7* 



(r) (r) 



5(r) |Ek 



W _ „(r) 



J'=l 



< n 



< 1. 



(A4) 



However, since =1 is integer, we con- 

clude that inequality (|A4[) is satisfied if and only if Eq. 
(|A2|) takes place. □ 



(A3) 
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